Impurity Entanglement Entropy and the Kondo Screening Cloud 
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The screening of an impurity spin by conduction electrons is associated with the formation of a 
large Kondo screening cloud, of size ^k- We study the quantum entanglement between a region 
of size r surrounding the impurity and the rest of the sample, (of total size TV) using Density 
Matrix Renormalization Group and analytic methods. The corresponding "impurity entanglement 
entropy", Simp, is shown to be a universal scaling function of r/^K and r/TZ. We calculate this 
universal function using Fermi liquid theory in the regime r. 
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Quantum entanglement of many body wave-functions 
is receiving much attention of late due to its quantum 
information (QI) theory applications 0, its importance 
to computational algorithms and the new perspective 
that it offers on critical behavior of condensed matter sys- 
tems 0, Q • Very recently, entanglement has been investi- 
gated in some quantum impurity models ■ The Kondo 
Hamiltonian 0, a famous quantum impurity model, was 
originally proposed to represent a single impurity spin 
interacting with conduction electrons in a metal. The 
Hamiltonian is: 



H 



dV[V/(-VV2m)i/; + JKS^{r)tl;\a/2)ij ■ S], (1) 



where ^{r) is the electron annihilation operator (with 
spin-index suppressed), the 5"° arc S = 1/2 spin oper- 
ators and the Kondo coupling, Jk > 0. Due to the 6- 
function in Eq. we may expand ipir) in spherical 

harmonics and only the s-wave couples to the impurity, 
yielding an effective one dimensional (ID) model. From 
spin-charge separation in one dimension, it can be seen 
that the universal low energy long distance physics of the 
Kondo model also occurs in the Ji — J2 family of gapless 
antiferromagnetic 5 = 1/2 Heisenberg spin chain models 
for J2 < J2 — 0.2412, with the coupling of the end sgin 
to its neighbors, J^, weaker than the other couplings 
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The Kondo coupling, Jj^ (or J'j^), is found to grow large 
at low energies under renormalization and the ground 
state is a highly non-trivial many body state in which 
the impurity spin is sometimes heuristically thought of 
as being screened by a conduction electron in a large 
orbital of size 



= v/Tk- 



(3) 



Here v is the velocity of the low energy excitations 
(fermions or spin-waves) and Tk is the energy scale at 




FIG. 1: (a) Schematic picture of the two regions A,B around 
a Kondo Impurity, (b) ID analogue of the Kondo problem, 
with TZ interacting qubits along an antiferromagnetic chain, 
(c) Total entanglement entropy, S{,J'if,r, TV) for a 102 site spin 
chain at JJ, with a J'^ = 0.3 Kondo impurity along with 
S{J'k — l,r — 1,7?. — 1) (no impurity). Uniform parts (solid 
lines) and the resulting Simp. 



which the rcnormalizcd Kondo coupling becomes large. 
For weak bare coupling, oc e"^"^^ where a is a con- 
stant. Many physical observables were shown to vary on 
the characteristic length scale, £,k [SIqI. Iio|. 

The purpose of this paper is to investigate the inter- 
play of the Kondo screening cloud and quantum entangle- 
ment occurring in the models Eqs. H1I2|) . In a QI context 
we could think of the three dimensional (313) fcrmionic 
model, Eq. (Q), as representing a single qubit coupled 
to a dissipative environment. Alternatively, the ID spin 
model, Eq. corresponds to a chain of qubits with the 
end qubit more weakly coupled than the rest. = 1/2 
chains are receiving considerable attention from the QI 
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community in the context of perfect state transfer 
In the 3D fermionic model, wc imagine separating space 
into a spherical region, A, of radius r surrounding the 
impurity spin and the region B, exterior to the sphere as 
in Fig. ^a). Tracing out the region B from the (pure) 
ground state density matrix, we obtain the reduced den- 
sity matrix p{r) for A. The entanglement entropy is de- 
fined as its Neumann entropy S = — tr[p(r) In p(r)] p^ . 
In the ID spin chain version, we consider the entangle- 
ment of the first r spins with the remainder of the chain 
as in Fig.^b). For the ID spin chain at J| wc show the 
resulting entanglement entropy in Fig. ^c) for 
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and with no impurity for J'j^ = 1. Note that, in both 
cases, there is a strong even-odd alternation of S with 
r. This effect was studied in 0| for the case of uniform 
coupHng, J'j, 



K 



1. In the present paper we focus on the 
uniform part, which can be defined at r 3> 1 by: 

5(r,7^) = Su{J'K,r,n) + {^ly Sa{J'k ,r,TZ), (4) 

where SuiJx,r,TZ) and 5'yi( J^, r, 7?.) arc slowly varying 
functions or r. There is also a strong alternation with 
the parity of TZ and we study both TZ even and odd. 
The uniform part consists of a bulk part, Suoir^TV) ~ 
(1 /6) ln[7^ sinfTrr /7^)1 -I- constant, independent of J'j^ 0, 
[13, [iJI and an impurity term that we define precisely as: 
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FIG. 2: Universal scaling plot of 5imp for fixed r/TZ, (a) for TZ 
even, (b) for TZ odd. DMRG results for the Ji — J2 chain at Jl 
for various couplings J^-. The lines marked 7r5_f(-/(12r) are the 
FLT prediction: Eq. 1101 . (c): the location of the maximum, 
(''/C/f )max, of Simp for odd TZ, plotted versus r/TZ. 



S,mp{JK.r,'JZ) = Su{J'K,r,TZ)~Su{l,r-l,n-l). (5) 

Here we have subtracted Su when the impurity is absent. 
See Fig. . There is also an impurity term in Sa but 
we do not address this here. Wc focus here on calculat- 
ing Simp for the spin chain model but the same universal 
scaling function also describes the additional entangle- 
ment entropy due to the impurity in the 3D Kondo model 
[l5| . We find that Simp{J'j^, r, TZ) appears to be a scaling 
function, depending only on the ratios of characteristic 
lengths when r, £_k ^ 1: 

StnipiJ'K, r, TZ) = Simp{r/£.K,r/TZ). (6) 

Note that there are actually two different scaling func- 
tions for the cases of TZ even or odd. We have verified this 
scaling form over a wide range of parameters using the 
Density Matrix Renormalization Group (DMRG) keep- 
ing from 256 to 1024 states (Fig. [2). In the hmit ^/^ < r 
we have calculated the universal scaling function explic- 
itly using conformal field theory (CFT) methods based 
on Nozieres' local Fermi liquid theory (FLT) 0,^3- 

Numerical Results. We have determined Simp for a 
range of values of J'j^ , r and TZ using DMRG on the spin 
chain model, Eq. In order to avoid well-known loga- 
rithmic corrections due to a marginally irrelevant (bulk) 
operator we tune J2 to the critical point J| « .2412, 
where the marginal coupling constant vanishes 0- We 
then extract Sjj from the DMRG data with high precision 



using a 7-point formula. In Fig. [21 we show a scaling plot 
of Sirnp for several fixed r/TZ allowing for a determination 
oi £,k{J'k) f-iid validating the scaling form, Eq. jnj. 

For odd TZ, S{0,r,TZ) - 5(1, r - l,7^ - 1) is exactly 
zero, since the decoupled impurity at j = 1 has no en- 
tanglement with the rest of the chain in the S ~ 1/2 
ground state in which it is fully polarized. This implies 
that Simp{0,r,TZ) = for 7?. odd. On the other hand, 
for TZ even, Simp{0, r, TZ) is not zero, as shown in Fig. 01 
Thus the presence of the impurity spin, even when there 
is no term in the Hamiltonian coupling it to the rest 
of the system, alters S. The reason this is possible is 
because, with TZ even and = 0, the ground state is 4- 
fold degenerate, corresponding to the impurity spin and 
the rest of the chain independently having 5^ = ±1/2. 
The ground state that we consider in this case is the 
spin singlet state. This is necessary to connect to the 
limit J'j^ 0+. S{0,r,TZ) is then clearly different from 
5(1, r — 1,7?. — 1) with a resulting non-zero Simp{0, r, TZ) 
for TZ even We expect that, in the limit r/TZ — > 0, 
Simpir / ^k) becomes the same for even and odd R, since 
adding one more site to an infinite chain should not have 
any effect on the entanglement entropy. Our data sug- 
gests that \im.r/^^^o\iToar/R^oSimpir/^K,r/TZ) = In 2, 
for either parity of TZ. We note that the presence of a 
term in the entanglement entropy at a conformally in- 
variant boundary fixed point, corresponding to the zero 
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temperature thermodyriamic impurity entropy, In g was 
pointed out in Ref. For the single channel Kondo 

model In g = In 2 at the weak coupling fixed point where 
the impurity is unscreened, and In g = at the strong 
coupling fixed point where it is screened. The thermo- 
dynamic impurity entropy decreases monotonically from 
In 2 to with decreasing T/Tk- The impurity entan- 
glement entropy behaves the same way with increasing 



stretching into region B. 
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FIG. 4: Typical valence bond configurations for TZ even (a) 
and TZ odd (b) where the unpaired spin is on the 7"^ site. 

FLT calculation. The non-interacting part of the ID 
(s-wavc) fcrmion Hamiltonian can be written as the in- 
tegral over r of spin and charge energy densities for left 
and right-movers, 'Hs,l/r and Hcl/r- At energy scales 
^ Tk the impurity spin is screened and does not appear 
in the low energy effective Hamiltonian. Irrelevant local 
interactions, at r = 0, not involving the impurity spin 
appear instead The most relevant of these can be 

written: 



(7) 



FIG. 3: DMRG results using spin-inversion for S'inip(0, r, 7?,) 
for TZ even. Inset: DMRG results for Simp (1, J", 7?,) for TZ 
odd and even. Data for TZ = 100, 101 and 7^ = 400, 401 
are indistinguishable. S'imp(l, 7?.) appears to vanish in the 
scaling limit r oo with r/TZ held fixed. 

A heuristic picture giving the general features of 
Simp{r/^K,r/TZ), for both even and odd TZ, can be de- 
rived from a resonating valence bond view of the ground 
state . We loosely identify Simp with In 2 x the prob- 
ability of an "impurity valence bond" (IVB) stretching 
from the impurity site into region B. (See Fig. 01) For 
J'j^ = and even TZ, all lengths for the IVB are equally 
probable so the naive picture predicts that Simp decreases 
linearly from In 2 to with increasing r/R, in rough 
agreement with the DMRG data in Fig. O We expect 
the typical length of the IVB to be (when < TZ), 
leading to the monotonic decrease of Simp with increasing 
r/^if , for 7^ even, and its vanishing as r/^x oo. For TZ 
odd, the ground state contains one unpaired spin. When 
^if = oo this unpaired spin is the impurity, implying no 
IVB and hence Simp = 0. As decreases the probabil- 
ity of having an IVB increases, since the unpaired spin 
becomes more likely to be at another site, due to Kondo 
screening, but the average length of the IVB, when it is 
present, decreases. These two effects trade off to give 
Simp a maximum when ~ R, in good agreement with 
Fig. El inset, at which point the probability of having an 
IVB becomes 0(1) but the decreasing average size of the 
IVB starts to significantly decrease the probability of it 



Various long distance, low energy properties can be cal- 
culated in an expansion in ^^/r and E^k/v at dis- 
tances r ^ ^/s" and energy scales E <^ v/^k- For in- 
stance, the zero temperature impurity susceptibility is 
^[T = 0) = CK/(4ti). In the limit < r, we can 
calculate Simp using the FLT interaction of Eq. (0 in 
lowest order perturbation theory in (^k- This calculation 
is easily done by CFT methods. We first express trp" in 
terms of the partition function, Zn{r) on an n-sheeted 
Riemann surface, TZn, with the sheets joined at the cut 
extending from r to 7?. 0, . The entanglement entropy 
is obtained by means of the replica trick: 

5 = -lim-^[Z„/Z"]. (8) 

(Z is the partition function on the ordinary complex 
plane.) The correction to Z„ of first order in is: 

/oo 
dT{ns,L{T,Q))n„. (9) 
-oo 

7is,L = T / (27r) where T{t, x) is the conventionally nor- 
malized energy- momentum tensor for the c = 1, free 
boson conformal field theory corresponding to the spin 
excitations of the original free fermion model. It obeys 
0: (T(r,0))TC„ = vl^ni'i.irf /[{vT + irf{vT - irf], with 
A„ = (c/24)(l - n"2). Both the integral in Eq. ^ and 
the limit in eq. ^ are trivial, yielding (see FiglSJ: 

S^mp = TT^/f /(12r). (10) 



Our FLT calculation of Simp is easily extended to finite 
Tt. In this case, (T) is given as above with the sub- 
stitution: VTzbir ^ (27^/7r) sinh[7r(uT ± ir)/(27?.)], and 
2r (27e/7r)sin[7rr/(7e)]. The integral in Eq. ® can 
again be done exactly, yielding: 
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[7^^/^/(127^)][l + 7r(l - r/7^) cot(7r^/7^)]. 



(11) 



We note that this result can be simply obtained from 
Suo ~ (1/6) ln[7^sin7rr/7?.)] by shifting both r and R by 
7r^A'/2 and Taylor expanding to first order in the small 
quantities ^if/r and ^k/ R- Excellent agreement between 
this result and DMRG data is shown in Fig.l^a) allowing 
for a determination of £,k{J'k)- 

In Figs. I2l5f a) we have determined £,k{J'k) by re- 
quiring good data collapse. It is an important consis- 
tency check on our scaling hypothesis to also determine 
^k{J'j) independently by more conventional methods. 
Using the beta- function to cubic order we find that p^ : 
oc exp(1.376/J^)/-y/ J^. [To derive the 1.376 factor 
we needed the spin- wave velocity, v k, 1.1699 and the 
coeflacient of proportionality between 5*2 + JIS's and the 
free fermion spin density at the impurity location, which 
was determined by comparing the end to end spin corre- 
lation function (decaying as \/Ti?) in the spin chain and 
\D free fermion models.] This formula is only valid at 
very weak coupling. We also derived Jl5| the Fermi liquid 
result: 
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(12) 

(where Cq, Ci and 62 are some constants) and used it to 
extract £,k{J'k) from DMRG results for the ground state 
energy. Finally, we compared the magnetization versus 
field of our spin model to exact Bethe ansatz results 
for the ID fermion model with length JZ, obtaining ex- 
cellent agreement and determining We show all 
our results for versus 1/J'x in Fig. I^b) observing a 
good agreement between the different estimates and, at 
small J'j^, with the weak coupling RG prediction. 

We note that, if we take region A to be the impurity 
spin only, then the resulting entanglement entropy can be 
expressed in terms of the ground state impurity magneti- 
zation, nii 



= -E±(l/2±: 



''imp 



) In (1/2 ± Wimp) 



(Here nii is only non-zero when the total number of 
spins is odd, in which case the ground state is a spin 
doublet.) rriimp was analyzed in [Sj for the ID fermion 
Kondo model and we find similar behavior for the spin 
chain model. This immediately determines the behavior 

of Simp HI- 

In conclusion we have defined an impurity entangle- 
ment entropy, calculated it using DMRG and CFT meth- 
ods and shown it to be a universal scaling function. We 
have also verified that the spin chain with a weak bound- 
ary link is indeed a realization of the Kondo model and 
determined accurately the corresponding Kondo length 
scale (Ck) as a function of this boundary coupling. 
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FIG. 5: (a) DMRG results for Si,-np{r)1l/ £,k for the spin chain 
at J2 for various values of J'^ in the Fermi liquid regime. 
The black curve is the FLT prediction, Eq. JTTJ. (b) The 
various estimates for compared with the RG result (^k oc 
exp(1.376/Jj^)/V^. 
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